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Abstract 

Food webs can be regarded as energy transporting networks in which the 
weight of each edge denotes the energy flux between two species. By inves- 
tigating 21 empirical weighted food webs as energy flow networks, we found 
several ubiquitous scaling behaviors. Two random variables Ai and Ci defined 
for each vertex i, representing the total flux (also called vertex intensity) and 
total indirect effect or energy store of i, were found to follow power law dis- 
tributions with the exponents a ~ 1.32 and /3 ~ 1.33, respectively. Another 
scaling behavior is the power law relationship, Ci ~ Aj, where t] ^ 1.02. This 
is known as the allometric scaling power law relationship because Ai can be 
treated as metabolism and Ci as the body mass of the sub-network rooted 
from the vertex z, according to the algorithm presented in this paper. Finally, 
a simple relationship among these power law exponents, rj = (a — l)/(/3 — 1), 
was mathematically derived and tested by the empirical food webs. 
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1. Introduction 

Scientists look for universal patterns of complex systems because su ch in- 



varia nt features may help to unveil the principles of system organization (jWaldrop 



19921 ). Complex network studies can not only provide a unique viewpoint 



of nature and society but also reveal ubiquitous patte rns, e.g., small world 



and scale free, characteristi c of va rious complex systems (j Watts and Strogatzl . 



19981 : [Albert and Barabasil . 120021 ). However, ecological studies have shown 
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that binary food webs, which depict trophic interactions in ecosy stems, refuse 

to become part of the sma ll world and scale free networks family ( iMontoya and Sole , 



2OO2I : iDunne et al.l . |2002| ). Although some common features, including "two 



degrees separation" , which means the ver y small average distance, are shared 
among food webs (IWilliams et al.l . |2002| ). other meaningful attributes such 



as degree distribution and clustering coe fficient change wit h the size and 



complexity (connectance) of the network (IDunne et al.l . |2002| ) 



Weight information of complex networks such as air traffic network or 
metabolism networks, etc., can reveal mo re unique patt e rns and features that 



are never found in binary relationships (IBarrat et al.l . |2004J : lAlmaas et al. 



2OO4J : iMontis et al.l . 120071 ). Food web weights have two different, yet c orre- 



lated , meanings in ecology . One is the strength of the trophic interaction (lEmmerson and Raffaelli . 



2OO4J : iBerlow et al.l . 120041 ): the other is the amount of energy flow. They are 



correlated because interaction strength is the per capita measure of energy 
flow. Interaction strength-based weighted food webs exhibit new features 
such as a relationship be t ween predator-prey body size ratio and in teraction 



strength flMcCann et all . 11998 



2OO2I : iBerlow et al.l . 12004 ). 



iWootton ■ 

Additionally, weights of food webs can also be denoted as the total 
amount of energy flow between two species when the whole system is in the 
stead y state . Stud i es of energy flow networks in ecosys t ems have a long his- 
torv (lOduml. Il988l: iFiiinl. Il976l: ISzvrmer and Ulaiiowiczl. 19871: 
Baird and Ulanowiczl . Il989l : iHigashi et al.l . Il993l : IPattenl . Il981 



Higashil . 11986 : 



19821). Many 



systematic in dicators were designed t o depict the macro -state of energy flow 



in ecosystem (iFath and Pattern Il999l : iFath et al.l . 1200 ll ) , of which some can 



not only reflect the direct energy flow s between species but al s o ind i rect ef 



fects and inter-dependence of species (IFath and Pattenl . Il999l : iFinnl . Il976l ). 
Although some important discoveries were made about the general struc- 
ture and function of ecological networks (iFath and B.C.I. Il9 98: Ha nnonl . ll973 : 



Levind . Il980l : iHannonl . Il986l : IPattenl . Il985l: luiano wiczl. Il986.. .19971). few focus 



on power law distributions and relations (lUlanowicz and Wolff. I . Il990l ). 

AUqmetric scaling is an important universal pattern of flow systems. 



Kleiberl (119321 ) found that the metabolism and body size of all species fol- 
low a ubiq u itous powe r law relationship, w ith an exponent around 3/4. 
West et al.l (Il997l ) and iBanavar et al.l ( 119991 ) explained this pattern as an 
emergent property of nutrient and energy transportation networks. This 
recognition encouraged people to realize that all qmetric scaling may be a 



universal feature for all transportation systems. iGarlaschelli et al.l (120031 ) 



extended Banavar's approach to binary food webs and found a similar al- 
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lometric scaling power law relationship. Although Garlaschelli's method as 
an algorithm had be e n app lied to various networks, including th e worldwide 



trade network (jDuanl, 120071 ) and tree of life ljHerrada et all l2008l ). it had sev- 



eral shortcomings. The first step of his algorithm is to obtain a spanning 
tree by cutting many edges in the original netw o rk so that a certain amount 



of information is lost (lGarlaschelli et al.l . |2003| ). lAllesina and Bodinil (j2005h 



improved this method by reducing the original network to a directed acyclic 
graph. Although less information is lost, cutting edges is still unavoidable. 
The second shortcoming of Garlaschelli's approach and Allesina's improve- 
ment is they can be applied to binary networks, but not weighted ones. 

This paper will combine the successful approaches in complex weighted 
networks and earlier studies on ecological flow networks to reveal the under- 
lying heterogeneities and universal scaling behaviors of food webs. The study 
is organized as follows. In section |2l the basic ideas and steps for obtaining 
Ai and Ci are introduced. Afterwards, we apply these tools on 21 empiri- 
cal food webs with energy flow information. Section 13.21 and 13.31 study the 
power law distributions of Ai and Cj. We extend Garlarschelli's approach 
to weighted food webs without cutting edges. The allometric scaling power 
law relationship between Ai and Ci is shown in section 13.41 A simple mathe- 
matical relationship among scaling exponents of power law distributions and 
power law relation is derived and tested using empirical food webs in section 
13.51 Finally, the ecological meaning of Ai and Ci, the distributions of flux 
matrix and fundamental matrix, consideration of node information, etc., are 
discussed in section HJ A simple example of our approach, a comparison to 
the existing methods, and the theorem regarding power law exponents are 
presented in the Appendix. 

2. Methods 

In this section, we outline the basic idea and mathematical definition of 
our method. One simple example showing how the approach works will be 



discussed in the Appendix A in detail. 



2.1. Flux Matrix 

An ecological energy flow network is a weighted directed graph that repre- 
sents relationships of ecological energy transfer. For a given graph, a matrix 
called flux matrix in this paper can be defined as representing the energy 
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flux between species. 



{N+2)x{N+2) 



:j}{N+2)x(N+2) 



(1) 



where fij is the energy flux from species i to j. Two special vertices represent 
the environment: vertex and vertex + 1. Vertex denotes the source of 
energy flow, whereas vertex + 1 represents the sink. We expect that the 
dissipative and exported energy will flow to vertex N + 1. Therefore, there 
are in total (A^ + 2) x [N + 2) entries in the flux matrix. 

2.2. Fundamental Matrix 

Suppose that the flow network is balanced, meaning that the total influx 
equals the efflux for each vertex i G [1,A^]. We can then define an. N x N 
matrix M from F follows, 



N+l 



m 



(2) 



k=l 



and the fundamental matrix can be derived as 



U = I + M + M"^ + ■ ■ ■ = ^M' = {I - MY 



(3) 



i=0 



where, / is the unity matrix. Any element Uij in U matrix denotes the 
influence i to j along all possible pathways . U rnatrix was first introduced 
in economic input-output analyses (iLeontiefl . Il95ll Il966l ) to indicat e the di 



rect a nd indirect effects of good flows in various ec onomic sectors. iHannon 



(I19731) was the f irst to apply this matrix to ecology( lFath and Patteru . Il999 ; 
Ulanowiczl - liooi ). 

Given the flux matrix and fundamental matrix, two vertex-related vari- 
ables, Ai and Cj, which will later be shown to follow power law distributions, 
are defined. 



2.3. Ai 

We can calculate the total flux through any given vertex i according to 
F. Th is value is also calle d node intensity in complex weighted network 
studies ( I Almaas et al.l . 12004 ). Because the network is balanced, we need only 
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calculate the efflux of each node as Ai, 



N+l 

A = $^/.i,V2G[l,iV] 



(4) 



2.1 Ci 

Another vertex-related index called Ci can be defined to reflect the total 
indirect effects or the total energy store of the sub-network rooted from vertex 



N N 

fc=i j=i 



/Uii)Uik 



(5) 



We will provide an explanation of the indicator Ci in Appendix A by a 
simple example. 

Ai is the total flow-through or intensity of vertex i. Ci is the total influ- 
ence of vertex i on all vertices in the whole network. Supp ose that of the 



19931 ). those passing 



many particles flowing in the network ( iHigashi et al. 
vertex i will be colored red. Ci would then be the total number of red parti- 
cles flowing in t he network. Actua l ly, the se two variables are extended from 
the approach of iGarlaschelli et al.l (|2003[ ) to calculate the allometric scaling 
of food webs 



see 



Appendix Bp 



2.5. Balancing the Network 

Sometimes the empirical network is not strictly balanced. 



algorithm, we can balance them artiflcially. Suppose X^yLo /i* X] 



To facilitate our 

N+l 



for 



=1 

vertex i. We can add an edge with the weight \flj\, flj = J2f=o fji~^^=i fij 
to connect the vertex i to A^-|- 1 or 0. If //^ > 0, the direction of this artiflcial 
edge is from i to N + 1. If f-j < 0, the direction is from to i. Normally, the 
artiflcial edges have very small weights because most empirical food webs are 
almost balanced already. 

2. 6. Power Laws 

After calculating the indicators of Ai and Ci for each vertex i, we will 
show that they follow the power law distributions in the high tails, which 
means that, 



P{Ai > x) ~ X 



l-a 



(6) 
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and 



Pia > ~ y'~^ (7) 
for given Ai and which are larger than given thresholds Xq, Uq (IClauset et al. 



20071 ). and where ~ represents "proportional to." The cumulative probability 



distribution curves will be shown and the power law exponents a, (3 calculated 
in the next section. 

Furthermore, we will show that Ai and Ci satisfy a power law relationship, 



a ~ Al 



(8) 



This relationship is also called the allometric scaling law because A^ rep- 
resents metabolism and Ci is the equivalent body mass of the sub-system 
rooted from vertex i 



see 



Appendix A). 



3. Results 

3.1. Dataset 

Twenty one food webs containing energy flow information from different 
habitats were studied(Table [T]). These food webs were obtained from an 
online databas^, and most are from published papers. In Table [H we list 
the name and the number of nodes (|iV|) and edges {\E\) in each web. The 
number of nodes does not include the "respiration" node, and the number of 
edges only counts the energy flows between species, and does not include the 
edges from (to) "input" and "output." The weights of edges in these food 
webs are energy flows whose values vary across a large range because the 
units and time scales of the measurements are very different. 

After studying the scaling laws of these food webs, we divided the results 
into two main sections. First, the power law distributions that reflect the 
heterogeneity of Ai and Ci are shown. Second, the allometric scaling power 
law relationship that depicts the self-similar structures of energy flows is 
discussed. 

3.2. Power Law Distributions of Ai 

We calculated the random variables Ai for each of the 21 empirical food 
webs. Four of them are selected to plot in Figure [TJ 



http://vlado.fmf.uni-lj.si/pub/networks/data/bio/foodweb/foodweb.htm 
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Table 1: Empirical Food Webs and Their Topological Properties 
(I AT I stands for the number of vertices of the network and is the number of 
edges. The webs are sorted by their number of edges.) 



Food web 


Abbre. 


N 


\E 


Florida Bay, Dry Season 


BayDry 


127 


1969 


Florida Bay, Wet Season 


BayWet 


127 


1938 


J- ivJiiLld D dy 


J- WjVWXOti 


197 


ic/OO 


Mangrove Estuary, Dry Season 


MangDry 


96 


1339 


Everglades Graminoid Marshes 


Everglades 


68 


793 


Everglades Graminoids, Dry Season 


GramDry 


68 


793 


Everglades Graminoids, Wet Season 


GramWet 


68 


793 


Cypress, Dry Season 


CypDry 


70 


554 


Cypress, Wet Season 


CypWet 


70 


545 


Mondego Estuary - Zostrea site 


Mondego 


45 


348 


St. Marks River (Florida) 


StMarks 


53 


270 


Lake Michigan 


Michigan 


38 


172 


Narragansett Bay 


Narragan 


34 


158 


Upper Chesapeake Bay in Summer 


ChesUpper 


36 


158 


Middle Chesapeake Bay in Summer 


ChesMiddle 


36 


149 


Chesapeake Bay Mesohaline Net 


Chesapeake 


38 


122 


Lower Chesapeake Bay in Summer 


ChesLower 


36 


115 


Crystal River Creek (Control) 


CrystalC 


23 


81 


Crystal River Creek (Delta Temp) 


CrystalD 


23 


60 


Charca de Maspalomas 


Maspalomas 


23 


55 


Rhode River Watershed - Water Budget 


Rhode 


19 


35 
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Figure 1: Cumulative probability of the variable Ai for four selected food webs with the 
best fitted lines. The fitted lines start from xq, and their slopes are listed in Table [2] 
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AiS of real food webs decay as power law ir i the high tail ; therefore, the 
fitted lines start from given lower bounds Xq ( IClauset et al.l . 120071 ) . These 
figures show that Ai follows the power law distribution. 

Th ere are obvious cutoffs in the tails that may be attributed to sampling 
effects ( iNewmanl . |2005| ). According to equation |6l the cumulative probability 
decays as x slowly. The probability of finding a higher Ai value (in the 
tail of the curve) is very small. Therefore, the number of samples in this 
interval becomes very few and statistical fluctuations are unavoidably large 
as a fraction of sample number . This phenomenon is obv ious in other fields 
such as income( 1Clementi et al.l . l2006l ). personal do nations (|Chen et al.l . l2009f ) 
and the number of species per genus of mammals ( iNewmanl . 120051 ). 

The scahng exponent a for ea ch food web was est imated according to 
the maximum likelihood approach (jClauset et al.l . l2007l ). The exponents and 
relative errors of power law fittings for all 21 food webs are listed in Table |2j 

In Table [21 xqs represent the lower bounds of the power law distributions. 
We normalized Xq by dividing the maximum A^ of each food web to avoid 
the large range variance of Xq among different webs because their units and 
the measurement time scales are very different. 



D is the statistic of the KS test dRousseau and Rousseaul . l2000l : lGoldstein et al. 



20041 ). Its value reflects the maximum distance between the cumulative prob- 



ability of real data and the fitted model. Therefore, the smaller D values 
indicate the better power law fitting, a is the quantile of the 95% confidence 
interval for different numbers of samples (INoetheii . Il967l ). and is only a ref- 
erence for p. If Z) i s sma ller than a, then we should accept the power- law 
hypothesis (INoetherl . Il967l ). From Table [2l we know that all food webs pass 
the KS test. In the last column, "No. of samples" means the number of 
sample points that are larger than xq and follow the power law distribution. 

By comparing different rows, we know that the food webs with more 
edges can be better described by power laws because their Ds are smaller. 
Further, the scaling exponent a and xo/xmin increase as the scale of the 
network decreases. All a values fall into the interval [1.23,1.77], with an 
average of 1.32. 

The power law distributions of AiS reflect the heterogeneities of energy 
flux. Few nodes possess high Ai values, while most nodes only share a small 
fraction of the energy flux. The exponent of power law reflects the degree 
of heterogeneity of the whole network. Therefore, larger food webs are more 
heterogeneous than smaller ones because their exponents are lower (Table 
[2]). Although the power law distribution of Ai can not give us concrete 
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Table 2: Parameters of the power law distributions for Ai 
(a is the power law distribution exponent of equation [H xq is the smallest value 
of Ai that follows power law, Xmax is the largest Ai; D is the KS statistic; a is 
the quantile of 95% confidence interval; and the number of samples is the total 
number of nodes following the power law distribution. The webs are sorted by 



their number of edges) 



Food web 




a 




max 


D 


cr 


No. of samples 


Baydry 


1, 


.25 


8, 


.45e- 


-006 


0, 


,09 


0, 


.13 


102 


Baywet 


1 


.23 


4, 


.24e- 


-006 


0, 


,08 


0, 


.13 


105 


Florida 


1 


.23 


4, 


.24e- 


-006 


0, 


,08 


0, 


.13 


105 


Mangdry 


1 


.24 


1, 


.55e- 


-006 


0, 


,06 


0, 


.16 


75 


Everglades 


1 


.23 


4, 


.20e- 


-006 


0, 


,12 


0, 


.20 


44 


Gramdry 


1 


.25 


9, 


.55e- 


-006 


0, 


,11 


0, 


.21 


41 


Gramwet 


1 


.23 


4, 


.20e- 


•006 


0, 


,12 


0, 


.20 


44 


Gyp dry 


1 


.28 


5, 


.95e- 


-005 


0, 


,10 


0, 


.20 


44 


Gypwet 


1 


.23 


2, 


,06e- 


•005 


0, 


,10 


0, 


.21 


42 


Mondego 


1 


.33 


3 


.87e- 


-004 


0, 


,11 


0, 


.26 


26 


StMarks 


1 


.35 


9, 


.36e- 


-004 


0, 


,17 


0, 


.21 


43 


Michigan 


1 


.32 


1, 


.80e- 


-003 


0, 


,18 


0, 


.31 


18 


Narragan 


1 


.26 


1, 


.lOe- 


-004 


0, 


,16 


0, 


.23 


33 


GhesUpper 


1 


.28 


9, 


.54e- 


-004 


0, 


,21 


0, 


.23 


32 


GhesMiddle 


1 


.29 


9, 


.47e- 


•004 


0, 


,17 


0, 


.24 


30 


Ghesapeake 


1 


.41 


8, 


.25e- 


-003 


0, 


,22 


0, 


.29 


21 


GhesLower 


1 


.77 


5, 


.79e- 


-002 


0, 


,20 


0, 


.33 


16 


GrystalG 


1 


.27 


1, 


.04e- 


-004 


0, 


,13 


0, 


.32 


17 


GrystalD 


1 


.25 


3 


■ lOe- 


-005 


0, 


,12 





.30 


19 


Maspalomas 


1 


.54 


2, 


.22e- 


-002 


0, 


,18 


0, 


.29 


20 


Rhode 


1 


.56 


1, 


.35e- 


-002 


0, 


,17 


0, 


.34 


15 
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Figure 2: Cumulative probability of the variable CiS for four selected food webs with the 
best fitted lines. The slopes of the fitted lines are listed in Table [3] 



information about each vertex (IFath and Pattenl . Il999l : iPattenl . Il98ll . Il982l ). 



it helps us to understand the network whole. 

3. 3. Power Law Distributions of Ci 

The same approach can be applied to CjS. The distributions of CjS for 
four selected food webs are shown in Figure O 

The curves of Ci distributions are very similar to the curves in Figure [H 
The estimated exponents and the KS test parameters are listed in Table [31 

Comparing Table [3] with Table |2l the exponents of Ci distributions are 
slightly higher than those of Ai distributions. The exponent /3 and KS test 
statistic D decrease with the scale of the network. The average exponent of 
these food webs is 1.33, with all Ci values falling into the interval [1.23, 1.84]. 
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Table 3: Parameters of the power law distributions for Ci 
{P is the power law distribution exponent of equation \7\ yo is the smallest value 
of Ci that follows the power law, Umax is the largest D is the KS statistic; a is 
the quantile of 95% confidence interval; the number of samples is the total 
number of nodes following the power law distribution. The webs are sorted by 

their number of edges) 



Food web (3 Vol Umax D a No. of samples 



Baydry 


1, 


.23 


2, 


.97e- 


-006 


0, 


,08 


0, 


.13 


107 


Baywet 


1, 


.23 


2, 


.71e- 


-006 


0, 


,08 


0, 


.14 


101 


Florida 


1, 


.23 


2, 


.71e- 


-006 


0, 


,08 


0, 


.14 


101 


Mangdry 


1, 


.25 


2, 


.03e- 


-006 


0, 


,07 


0, 


.16 


68 


Everglades 


1, 


.23 


2, 


.16e- 


-006 


0, 


,10 


0, 


.20 


44 


Gramdry 


1, 


.25 


4, 


.57e- 


-006 


0, 


,10 


0, 


.21 


42 


Gramwet 


1, 


.23 


2, 


.16e- 


•006 


0, 


,10 


0, 


.20 


44 


Gyp dry 


1, 


.28 


4, 


.88e- 


-005 


0, 


,10 


0, 


.20 


44 


Gypwet 


1, 


.25 


3, 


,13e- 


•005 


0, 


,11 


0, 


.21 


39 


Mondego 


1, 


.29 


9, 


.61e- 


-005 


0, 


,12 


0, 


.24 


30 


StMarks 


1, 


.42 


1, 


.27e- 


-003 


0, 


,16 


0, 


.21 


40 


Michigan 


1, 


.30 


7, 


.22e- 


-004 


0, 


,17 


0, 


.31 


18 


Narragan 


1, 


.24 


1, 


.22e- 


-004 


0, 


,17 


0, 


.23 


33 


GhesUpper 


1, 


.28 


6, 


.03e- 


-004 


0, 


,20 


0, 


.23 


32 


GhesMiddle 


1, 


.26 


3 


.88e- 


•004 


0, 


,19 


0, 


.24 


30 


Gliesapeake 


1, 


.61 


2, 


.41e- 


-002 


0, 


,18 


0, 


.33 


16 


GhesLower 


1, 


.84 


5, 


.02e- 


-002 


0, 


,22 


0, 


.35 


14 


GrystalG 


1, 


.28 


1, 


.04e- 


-004 


0, 


,13 


0, 


.32 


17 


GrystalD 


1, 


.23 


1, 


.35e- 


-005 


0, 


,12 





.29 


20 


Maspalomas 


1, 


.56 


1, 


.63e- 


-002 


0, 


,19 


0, 


.29 


20 


Rhode 


L 


.52 


8, 


.93e- 


-003 


0, 


,23 


0, 


.32 


17 
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We have studied the heterogeneity of CjS node by node. The nodes with 
high Ai values always have high Cj values, indicating a possible positive 
correlation between Ai and Cj. 

3.4- Allometric Scaling Relations 

The similarity between Figure H] and Figure [T] shows that there must 
be some connections between Ai and Cj. Allometric scaling of these flow 
networks revealed that the relationship between Ai and Ci is actually a power 
law. 

As shown in Figure [3l the sample points aggregate around their fitted 
lines very well. This relationship is ubiquitous for all 21 food webs as shown 
in Table H 

We used the minimum square error method to find the best-fitted line 
(Table H]). R'^s were larger than 0.9 for all food webs except CrystalC, Crys- 
talD and Rhode, whose scales are very small (|A^| < 23). The R^s and 
exponents decrease with the scale of the network because the statistical sig- 
nificance decreases as the number of samples declines. All exponents t] fall 
into the interval [0.83, 1.05]. The mean value of 77s for these food webs, except 
CrystalC, CrystalD and Rhode, is 1.02. 

We also show the Ai and Ci values of root nodes for all food webs, and 
fit them with a line on the log-log plot (Figure H]). These power law relations 
reflect the self-similar nature of the weighted food webs. 

3. 5. Relationship of Scaling Exponents 

As we have shown, Ai and Ci are random variables following power law 
distributions with scaling exponents a and (3, respectively. They also follow 
a power law relationship with the scaling exponent r/. Is there any universal 
relationship among a, /3 and 77? 

Actually, for any random variables following power law distributions and 



power law relations, we can prove a mathematical theorem (see Appendix C ). 
According to this theorem, the power law exponents a, /3, rj of the food webs 
should also satisfy equation IC.ll We tested this hypothesis by calculating 
{a — l)/{f3 — 1) — rj for each of the 21 empirical food webs to obtain Figure 

El 

From Figure El we know that most exponents of food webs satisfy this 
relation, and that the errors become larger as the scale of the network de- 
creases. As the scaling behaviors we are studying are statistical properties. 
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0.001 0.1 10 1000 10-= 0.001 0.1 10 1000 

Ai A,- 
Florida Mangdry 



Figure 3: AUomctric Scaling relationship between AiS and CiS of four selected food webs 
with the best fitted lines 
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Table 4: Allometric scaling of Empirical Food Webs 
(The second column lists rjs of each food web with the errors; The webs are 
sorted by their number of edges) 



Food web 




7] 






i?2 


Bay dry 


1 


.01±0, 


.01 


0, 


.9946 


Baywet 


1 


.02±0, 


.01 


0, 


.9946 


Florida 


1 


.02±0, 


.01 


0, 


.9946 


Mangdry 


1 


.01±0, 


.01 


0, 


.9967 


Everglades 


1 


.02±0, 


.01 


0, 


.9992 


Gramdry 
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Figure 4: Ai, d plot for root nodes of 21 food webs. The slope (the power law exponent) 
is 1.02±0.02 

the significance of tliese regularities will increase with the number of sam- 
ples. Therefore, scaling behaviors are more obvious and accurate for large 
scale networks because larger webs have more sample points. 

4. Discussion 

4.I. Ecological Meaning of Power Law Exponents 

As demonstrated above, food webs as energy transportation networks 
always follow power law distributions and relations. Three important ex- 
ponents (a,/3 and t]) are derived from these power law regularities. The 
question of whether these exponents carry ecological meaning naturally fol- 
lows, and at first, the three exponents all reflect integral properties of whole 
networks, a describes the heterogeneity of first passage energy flux distri- 
butions among vertices. The heterogeneity decreases with a. Therefore the 
distributions of energy flows are more uneven in large food webs than the 
small ones. From table |2l we also know that all a values fall into the interval 
[1.23, 1.77]. According to the features of power law distributions, the means 
and variances of power law random variables with exponents smaller than 
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Figure 5: The differences between t] and (a — l)/(/3 — 1) are calculated for all food webs 
to test the exponents relationship (equation IC.ip . The food webs are sorted according to 
their number of edges 
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2 are divergent (jNewmanl . |2005| ). Therefore, energy flux on food webs has 
no characteristic value. It is meaningless to find a specif ic species with the 
average energy flux as the representation of other species (jNewmanl . |2005| ). 
The allometric scaling r elatio n describes the self-similarity of flow net- 



works. iGarlaschelli et al.l (|2003[ ): lAUesina and Bodinil ( 120051 ) pointed out 



that allometric scaling exponents describe the transportation efficiency of 
binary food webs because Cj is treated as the cost of transportation. The 
range of these exponents is between 1 (most inefficient network) and 2 (most 
efficient network). However, we believe that the exponent t] discussed in this 
paper does not describe the efficiency of the whole network. As pointed out 



in section [2] and Appendix A , Ci can be understood as the energy store by 
the system, rooted from i but not the cost of the transportation. Thus, the 
food web with higher rj can store more energy with the same consumption 
of metabolites (Ai). Therefore, we believe that the food webs with higher 
7] are more capable of storing energy by means of cycling the flows in the 
network. In Table HI we see that the networks with larger scales have larger 
rj values. Consequently, food webs can increase their ability to store energy 
by increasing their complexity. 

Further, the range o f exponents r? is not s i mply [1 , 2] (see Table HI). As 
Garlaschelh et al.l (l2003h : lAUesina and Bodinil (l2005l ): iBanavar et "allll999h 
pointed out, the range [1,2] is only suitable for spanning trees or directed 
acyclic graphs of the original binary food webs. However, our method con- 
siders more ingredients, including the energy flux as the weight of edges, the 
loop structures of energy flows, and the heterogenous energy dissipation of 
each node, than the mere topology of the food webs with homogenous nodes. 
That is the reason why the exponents are out of the range [1, 2]. 

The exponent /3 also describes the heterogeneity of indirect effects. It is 
determined by exponents a and /3 via the theorem mentioned in section 13.51 



4-2. Flux Matrices and Fundamental Matrices 

As discussed in section [21 Ai and Ci are defined according to the fiux 
matrix and fundamental matrix. Therefore, the scaling behaviors of the food 
webs are determined by the matrices. The properties of these matrices may 
help us to understand the origin of the scaling behaviors. 

The elements in flux matrices also follow power law distributions, with 
an average exponent 1.46 (see Table |5]). We hypothesis that this power law 
determines the power law distribution of Ai. 
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Table 5: Parameters of the power law distributions for fij 
(a is the power law distribution exponent; Fq is the smallest value of fij that 
follows the power law; Fmax is the largest fij; D is the KS statistic; a is the 
quantile of 95% confidence interval; number of samples is the total number of 
edges following the power law distribution. The webs are sorted by their number 
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Fundamental Matrix Fundamental Matrix 

Figure 6: Cumulative distribution of fundamental matrices of four selected food webs 

However, unlike other variables, the distributions of fundamental matrices 
are not power laws but rather more like log-normal because the tails of the 
curves decline quickly (see Figure |6]). We also studied all of the fundamental 
matrices of 21 empirical food webs, and noted that very few could pass the 
KS test. 

We presume that the calculation of the fundamental matrix in equation [3] 
needs infinite operations on F matrix (see section [2]). As a result, the noise in 
F matrices is enlarged and accumulated in the tails of the distribution curves 
of U matrix. However, the means by which the non-power law distribution 
of fundamental matrices determines power law distributions of Ci and the 
allometric scaling power law relationship is an interesting problem for future 
studies. 
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4-3. Information on Nodes 

One of the weak points of our study on allometric scaling power law re- 
lationships is that the exponents are very close to 1. However, in this case, 
"allometry" just means the non-linear relationship between two variables, 
so, the relationship between Ai and Ci cannot be rigorously defined as an 
allometric scaling relationship. Because the calculations of fundamental ma- 
trices and CiS are always based on linear algebra, the results are close to 
linear relationship. One possible way to mend this weak point is to further 
consider the information available about nodes. 

Indeed, much information on nodes, i.e., each species in the web, is ig- 
nored in this work. The biomass as the weight of each node is available for 
many food webs. According to the definition, Ci is simply the energy store 
of the sub-system rooted from the vertex i. Therefore, biomass information 
should be included in Ci because a large part of energy will flow into the 
species node stored as biomass. It is possible that a new approach of cal- 
culating Ci including the biomass information for all species may break the 
linear relationship between Ai and Ci. 

Another important node characteristic is the body size of a given species. 
The metabolic theory predicts that species body size of the species c an no t 



only determine metabolism, life span, and birth rate, etc. (IBrownl. l2004f) . 



but a lso play an important role in energy flows and food webs ( 1 Cohen et al. 



20031 ) . An integrated theory of weighted food webs based on energy flow 



networks should contain body size data. 



5. Concluding Remarks 

This paper presents a new approach to reveal the scaling natures of 
weighted food webs as energy flow networks based on flux and fundamental 
matrices. The Ai,Ci distributions and the relationship between them always 
follow power laws. The power law exponents a, (3 and t] satisfy a relation- 
ship, = (a — — 1) as proved by the theorem [TJ Power law exponents 
consistently change with network scales. 

We note that the allometric scaling exponent does not reflect the trans- 
portation efficiency of networks but rather the capability of storing energy, 
which is very different from previous studies. We also investigated the dis- 
tributions of flux matrices and fundamental matrices, and suggested that 
biomass information should be incorporated into future studies. 
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Appendix A. A Simple Example 

To understand the method introduced in the section [21 let's look at a 
simple example (Figure rA.7p . 




Figure A. 7: An example to illustrate our method on deriving AiS and C^s 
(a) is the original flow network; (b) is the balanced network from (a), the dashed 
arrows are artificial edges pointing to the sink node; (c) is the Markov chain 
calculated according to (b). 



As shown in Figure \A.7\ the balanced flow network and the derived matrix 
M can be obtained step by step according to the method described in section 
m The fundamental matrix U can then be calculated for this simple network. 



U = I + M + M^ + ■■ 



( 1 






3/5 
42/41 
12/41 
3/164 

3/82 



7/20 
7/82 
42/41 
21/328 
21/164 



1/5 
14/41 
4/41 
165/164 
1/82 



2/5 \ 
28/41 

8/41 
21/41 
42/41 ) 



(A.l) 



Any entry mjj in the matrix M is rnerely the probability of one particle 
flowing from vertex i to vertex ?'( lBarberl . 119781 ). Furthermore, any entry (z, j) 
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in M ■ M represents the probabihty of a particle flowing from i to j along 
any path in 2 steps. M ■ M ■ M represents the probabilities after 3 steps, etc. 



(iHigashi et al.l . Il993l ). Thus, the matrix U simply takes in consideration all 



transfers of particles along all possible paths. 

Now, we will show how to calculate Ai and Ci for vertex 2. According to 
the definition, Ai is the total flow-through of vertex i, so A2 = Yl^j=i /a.j = 60. 

Suppose many particles are flowing in the network. They will be colored 
red once they flow through vertex 2. These particles will keep their color and 
flow around the whole network along all possible pathways. Hence, the total 
number of red particles in the network is just C2, which is computed as. 



C2 = G2J2 «2fc = ((100 X 3/5 + 0)/(42/41)) ^ U2k = 125 (A.2) 



k=l 



k=l 



where the term G2 = Yl^j=i fojUj2/u22 in equation IA.2I is the total num- 
ber of new particles which are colored red by vertex 2 in each time step. 
Ylf=i fojUj2 is the number of particles that flow into the system from the 
environment to the vertex i along all possible pathways at each time step, 
with foj = (100, 0, 0, 0, 0) in this example. By dividii ig by the term Uu t o 
derive Gi one avoids double counting the red particles (lHigashi et al.l . Il993l ). 
Thus, Ci = Gi Ylk=i "^ik is the total number of particles that have been col- 
ored red and flow to other nodes along all possible pathways at each time 
step. 

If we treat the red particles flowing in the network as a metabolic sub- 
system, we can calculate its all ometric scaling relations hip as Garlaschelli 
has done for binary food webs (j Garlaschelli et al.l . 120031 ). Thus, Ai is the 
metabolism and Ci is the energy store or body mass of the sub-system. 
Indeed, Garlaschelli's approach can be recovered by our method, as shown 
in the next section. 



Appendix B. Comparisons to Existing Approaches 



In this se ction, we will compar e our a pproach to lGarlaschelli et al.l ( l2003l )'s 
method and lAllesina and Bodinil ( l2005l )'s method. 

Figure IB.8( a).(b) shows how Garlaschelli's approach can be applied to 
a hypothetical food web to calculate Ai and Ci for each vertex. At first, 
a spanning tree is constructed from the original food web (Figure IB. 81 (a)) 
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Figure B.8: Calculations of allometric scaling of a hypothetical food web by various 
methods. 

(a) is a hypothetical food web (The letter in each vertex is its index). The black 
vertex is the root; (b) is a spanning tree of the original network (a). Ai and 
Ci are denoted inside and beside vertex i; (c) is the implicated flow network of 

(b) , the numbers beside edges are flux. The dashed lines are additional edges; 
(d) is a directed acyclic graph of the original network, the numbers are AiS and 
CjS calculated by the method of (lAllesina and Bodinil . l2005l ): (e) is a constructed 
flow network according to (d); (f) is the network with numbers of AiS and CjS 
calculated by our method according to the flow structure of (e). 
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by cutting edges. That way, each sub-tree rooted from any vertex can be 
viewed as a sub-system of the spanning tree. For example, the sub-tree with 
three vertices (b,f,g) rooted from the vertex b is a sub-system of the spanning 
tree. Ai is the total number of vertices involved in this sub-tree and Ci is the 
summation of AiS for each vertex in this sub-tree. Therefore, in this example. 
Ah is 3 and Cf, is 6. Finally, the universal allometric scaling relationship of AiS 
an d C,;S, with an exp o nent around 1.3, was found for all food webs, according 



to 



GarlascheUi et all (120031). 



Garlaschelli's method was inspired by iBanavar et al.l ( 1l999l )'s model to 
explain the Kleiber's law (See Figure iKST c)). The spanning tree is simply 
Banavar's optimal transportation network. Thus, energy flows into the whole 
system from the root along the links of the network to all nodes. Suppose 
that each node would consume 1 unit of energy in each time step. A finx 
with 1 unit representing the energy consumption by each node should then be 
added to the original spanning tree. In Figure iBlST c). the energy dissipation 
by each node is added as a dotted line. As a result, A^ of each node is just 
the total influx of this node. Ci is the total flux (the total number of red 
particles colored by i) of the sub-tree rooted from i. Essentially, calculation 
of allometric scalings using Garlaschelli's approach is based on this weighted 
flow network model. Therefore, our algorithm can derive the exact same 
val ues of Aj and Cj for t he flo w network (Figure iRST c)). 



AUesina and Bodinil ( 120051 ) extended Garlaschelli's method. First, the 



original network is converted to a directed acyclic graph (DAG) as shown 
in Figure iKST d). then the sub-network originated from vertex i is identified 
as the set of vertices that have at least one path from i. Therefore, vertices 
b,c,d,f,g,h belong to the sub-network rooted from b because they are all 
connected with b. A^ is the number of vertices in the sub-network, and Cj 
is the summation of all AiS in this sub- network as shown in Figure [K8] (d). 
Because AUesina and Bodini's approach is not based on weighted flows, it 
cannot be covered by our approach. 

However, we can construct a balanced flow network according to the 
original network as shown in Figure IB. 8( e). The information of edges is 
added. Our approach can be applied to this flow network to calculate Ai 
and Ci values (Figure iRST f)). Comparing Figure iRST d) to (f), we find their 
AiS are almost the same, except vertices a and h. To balance the network, 
more flows are added in vertex h and a, so their Ai values are larger than 
those in Figure iKST d). This modification can not only influence Ai values 
but also CjS. The vertices in networks IB. 8( d) and (f) therefore have different 
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Ci values. 



Although our approach requires weight information, it can be extended 
to more general flow networks, even those with cycles and loops. Also, as 
demonstrated by our approach, Ci merely means the energy stored in the sub- 
system rooted from vertex z, which provides much clearer and more significant 
ecological meaning than previous works. 

Appendix C. A Theorem about Power Law Exponents 

Theorem 1. Suppose X and Y are two random variables following power 
law distributions. Their density functions, p{x) ~ and p{y) ~ y~l^ hold 
for any positive x > Xo,y > yo, where Xq and yo are lower bounds of X and 
Y . Additionally, X and Y satisfy a power law relation, Y ~ X^ , then the 
exponents C(,f3,rj have following relationship: 



Proof. Because X and Y follow power law distributions, p{x) = cx~", where 
c is a constant that satisfies the normalization condition. And Y = kX^, 
where /c is a constant, so, for any y > kxQ, 



P{Y >y} = PikX"^ >y} = P{X > {y/ky/"} = / cx-"dx (C.2) 



Let t = kx"^, so X = {t/ky/'i, dx = {l/ri)k-^/H^^/'^^~^dt. Take it into equation 



P{Y>y}= / (cA;("-^)/Vr7)t(i-")/''-Mt = / k't^^-'^^'^^-^dt (C.3) 



where k' = ck^"" i)/'?//^ is a constant. Because Y follows power law distribu- 
tion. 



= (a - - 1) 



(C.l) 







y 



(C.4) 



Compare equation IC.3I with equation IC.4t we know that. 



/3 = {l-a)/r]-l 



(C.5) 
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Finally, 



= (a - - 1) 



(C.6) 
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